Kinks in the periodic Anderson model 
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Recently, dynamical mean field theory calculations have shown that kinks emerge in the real part 
of the self energy of strongly correlated metals close to the Fermi level. This gives rise to a similar 
behavior in the quasi-particle dispersion relation as well as in the electronic specific heat. Since 
f-electron systems are even more strongly correlated than the -hitherto studied- d-electron systems 
we apply the dynamical mean field approach with the numerical renormalization group method as 
impurity solver to study whether there are kinks in the periodic Anderson model. 

PACS numbers: 71.27. +a, 71.10.Fd 



In a general sense, kinks are an abrupt but in reality of- 
ten smooth crossover between two parameter regimes of 
a given physical system. Kinks are well known to result 
from the interaction between fermionic degrees of free- 
dom and external bosonic modest The coupling to collec- 
tive excitations such as phonons influences the electronic 
dispersion. This results in kinks inside an energy range 
of the order of the Debye frequency 10 p, centered around 
the Fermi level Ep. Typically these kinks are found at 
±40-70meV£ 

Due to the restriction to this rather small energy-range 
determined by up, the coupling to phonons cannot be the 
source of the observed high-energy kinks in the disper- 
sion at energies > 80mcV. Such high-energy kinks have 
been found, among others, in cuprates^ In these cases, 
a mechanism not depending on interactions with exter- 
nal bosonic degrees of freedom has to be the microscopic 
origin. 

It was recently discovered^ that such high-energy 
kinks emerge as an intrinsic feature of strongly corre- 
lated metals, in the real part of the self energy. The 
mechanism^ indeed requires no additional coupling to ex- 
ternal collective excitations and the corresponding kinks 
can arise at energies as high as a few hundred meV. 

A mathematical understanding^ can be gained on the 
basis of the Hubbard model within dynamical mean- 
field theory (DMFT)£"£ In the correlated metallic regime 
with a characteristic three peak spectral function, the 
kink energies can be shown to be dependent only on the 
renormalization strength Zpp and the non- interacting 
bandstructure£ For a Bethe-lattice with bandwidth W = 
2D the kinks are located at 



= Z FL (V2-1)D. 
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The effective dispersion before and after this kink follows 
directly by a renormalization of the free dispersion e k and 
is given by 



c_ + Zcp(e k - no) 
c + + Z£ p (e k - no) 



for uj < —uj*_ 

for — UJ*_ < LU < LO* + 

for ui > lu*_ 



Here, [1q denotes the chemical potential in the non- 
interacting case and Zqp is a second, weaker renormal- 
ization factor for quasiparticlc energies beyond u>*. 

The kinks in the real part of the self energy were also 
shown to reflect in a maximum in the susceptibility^ and 
to result in corresponding kinks in the low temperature 
electronic specific heat Cv{T)^ For the latter, one can 
estimate a kink temperature 
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Unfortunately, this kink temperature is for transition 
metals usually very large, i.e. of O(10QQK). At such high 
temperatures the specific heat is dominated by its phonon 
contribution, making an analysis virtually impossible. 
However, there is the important exception of L1V2O4, the 
first heavy Fermion system with d-electrons^ and a kink 
temperature of about 10K, confirming the theory^ Nat- 
urally, one would hence look at f-electron systems with 
a similarly low energy scale. However, at present it is 
unclear whether the electronic kinks of Ref. [H are to be 
expected at all for f-electron systems. They do exist for 
the Hubbard model^ but not for a single impurity Ander- 
son model with a constant conduction electron density of 
states^ Hence we ask ourselves: Are there purely elec- 
tronic kinks in the periodic Anderson model, the most 
fundamental model for f-electron systems? 

This paper addresses directly this question and shows 
the emergence of kinks in the real part of the self energy 
ReS(w) of the periodic Anderson model and the resulting 
effective energy- momentum dispersion E k . The outline 
is as follows: In Section U we first introduce the peri- 
odic Anderson model used for the analysis as well as the 
DMFT calculations themselves, employing the numerical 
renormalization group (NRG) as an impurity solver In 
Section |H] the results for the self energy £(cj), dispersion 
.Ek, and specific heat CV(T) are discussed. The main re- 
sults are summarized in Section IIIII 
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FIG. 1. (Color Online) Depiction of the PAM. The interacting 
f- level (ef,U) of each site couples via the hybridization V to 
the conduction band £k- Direct hopping between different f- 
levels is not possible. 



I. MODEL AND METHODS 

The focus of this work lies on systems of strongly 
correlated f-electrons. These are confined to very nar- 
row orbitals and interact with a local Coulomb repulsion 
U . Together with non-interacting conduction electrons 
these are the ingredients of the periodic Anderson model 
(PAM, Figure [T]). The corresponding Hamiltonian reads 



Hp am = e k a L a ka + £ fJ2 A>/i<r 

+ £K «/ kCT + ikVkJ (4) 

i ia 

Here, the operators a, a + represent a non-interacting 
conduction band with dispersion £k that can be thought 
of as a combination of s-, p- and d-bands, whereas the op- 
erators /, / + stand for the localized f-electrons with con- 
stant non-interacting energy Sf. Each site of this model 
consists of a conduction band site as well as a f-orbital, 
which are hybridized with each each other by strength 
V k . A detailed study of the PAM with DMFT(NRG) can 
be found, e.g., in Ref. [l4|. For the rest of this paper, the 
chemical potential fi is set to zero for the sake of simplic- 
ity. 

As for the non-interacting case, the Hamiltonian Eq. 
(U) can be more conveniently written as a 2 x 2 or- 
bital matrix for each k-point (in terms of conduction 
and f-electron orbital). Due to the interaction U, the 
f- level is modified by a self energy which is 

momentum- independent within DMFT. The correspond- 
ing one-particle Green's functions reads 



G k (w) = (col - ffk)- 1 = 



e f - EH 



W - £k 



(5) 

For the calculations carried out in this work an energy 
e c = -j^j- ^ k £ k is considered in Eq. (UJ) ^ which can be 
interpreted as the center of mass of the conduction band 
and acts as an energy-shift of the conduction electrons. 
Additionally, for computational reasons, the momentum 
dependency of the hybridization is neglected for the rest 




FIG. 2. (Color Online) Upper panel: local f-electron spec- 
tral functions of the PAM on a simple cubic lattice for 
ej = -0.5, U = 1, V 2 = 0.1, W = 2 and different values 
of e c = 0.1,0.3,0.5 (green/light gray, red/dark gray, black). 
Lower panel: corresponding real part of the self energy. In- 
sets: zoom in around the Fermi level at lo — 0. The kinks are 
indicated by arrows for the case e c = 0.5. 



of this work Vk V. This leads to the following expres- 
sions for the local Green's functions 



G/M = 



d 3 k 



1 



(2tt) 3 CM - £k 
1 



(6) 



V 



with £(u>) = uj 



V 2 



(7) 
(8) 



Eq. (J7J will serve as self-consistency relation for the 
DMFT calculations in this work. 



DMFT self-consistency cycle 

The complex many-body problem is approximated by 
the DMFT self- consistency cycled For the PAM imple- 
mentation, we need to solve the integral in Eq. (|6|). For 
the Bethe lattice with semi-elliptical density of states 
N Bethe (e) 



analytically, yielding 



£ , the integral can be calculated 



G c (lj) 
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FIG. 3. (Color Online) Upper panel: local f-electron spectral 
function for U = —2sf = 1, e c = 0.5, V 2 = 0.1 and different 
bandwidths W = 2,4,6 (black, red/dark gray, green/light 
gray) for a simple cubic lattice. Lower panel: corresponding 
real part of the self energy. Insets: larger frequency window. 
The kinks are indicated by arrows for the case W = 2. 



For the simple-cubic lattice, integral Eq. ([5]) can be 
rewritten in the formic 



G c {u) = -^P(z), (10) 

with the abbreviation z(w) = — gCTZJ) (W ~ is the 
bandwidth). Here, the function P(z) is equivalent to 
a product of two hypcrgcomctric functions of the type 
2 Fi(l/2, 1/2; l;fc 2 ), namely 



P(z) 



The abbreviations k±,xi and £2 arc defined by 



(11) 



k% = -±-x 2 V4^-^(2-x 2 )Vl zr ^2 (12) 
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Unfortunately, to treat the PAM with DMFT(NRG) 
involves overcoming numerical dimcultic o 13 ' 14 which af- 
fect the stability of the calculations (the NRG step in 
particular) as well as the convergence behavior of the 
DMFT-loop. Specifically, to arrive at fixed and stable so- 
lutions of the PAM for the considered parameter regime, 
it was necessary to introduce a small imaginary shift of 




FIG. 4. (Color Online) Real part of the self energy for the 
Bethe lattice with U = -2s f = 1, e c = 0.5, W = 2 and 
various V. Right inset: magnification around the Fermi en- 
ergy. Left inset: frequency of the kink vs. inverse hybridization 
strength. 



the real frequencies 
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i<5, which takes care of 



possible delta peaks in the Green's function Eq. (J7J or 
hybridization function, and to make use of Broyden's 
method of convergence stabilization . 17 ' 18 The logarith- 
mic NRG discretization has been taken to be A = 2 with 
checks for A = 1.8, 1.9, 2.0 yielding very similar results. 



II. RESULTS 

It is known that due to the hybridization of conduc- 
tion and f-electrons, the PAM at half filling (sf = —U/2 
and £ c = 0) represents a Kondo-insulator in the param- 
agnetic phase studied throughout this paper However, 
searching for kinks, a metallic configuration is necessary, 
and one with a well defined three peak spectral func- 
tion desirable. These requirements can be achieved by 
breaking the particle-hole symmetry in such a way that 
the energy level of the conduction electrons is increased 
while the f-level is held symmetric , 13 ^ see Fig. [2] (upper 
panel). 

The effect of the conduction band shift e c on the real 
part of the self energy is illustrated in the bottom panel of 
Fig. [2] For w < 0, the overall behavior of the self energy 
is the same as one would expect for a strongly correlated 
metal. The real part shows a basically linear behavior 
for small energies and eventually reaches its maximum. 
After that, it falls off rapidly and ultimately converges to 
a constant. 

For u> > 0, on the other hand, S(w) experiences 
the consequences of the hybridization. The hybridization 
gap, which moves to higher frequencies for increasing e c , 
is reflected in the imaginary part of the self energy as a 
second minimum of |Im£(u;)| (besides oj = 0; not shown). 
In the Kramers-Kronig related Re£(u;) it shows up as an 
inflection point in Fig. [2J 
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log(A f (co,k)+A c (co,k)) log(A f (<D,k)+A c (<o,k)) 




(a) (b) (c) 

FIG. 5. (Color Online) (a) Overall spectral density (A/(k, w) + A c (k, ui)) on a logarithmic scale, for a simple cubic lattice with 
U = -2s j = 2, s c = 0.5, V 2 = 0.2, and W = 3.5. (b) magnification around the Fermi energy. The curves along the maxima 
of (a) and (b) represent the dispersion E^. In (c) this dispersion (solid red curve) is compared to the non-interacting one 
multiplied with renormalization factors Zcp = 0.202 (green dotted) and Zfl = 0.0689 (blue dotted) extracted from a fit of 
the corresponding self energy. The dotted curves fit the DMFT(NRG) result very well, hence indicating kinks in the dispersion 
of the PAM. 



Compared to these coarse features, kinks are fine struc- 
tures which have been overlooked in the Hubbard model 
prior to Refs. i and H, and in the periodic Anderson 
model up to this point. A closer inspection of the self- 
energy shows a kink for u> < close to the Fermi en- 
ergy, see the inset of Fig. [3J particularly well visible 
for e c = 0.5. Upon increasing the bandwidth W and 
hence decreasing the density of states of the conduc- 
tion electrons at the Fermi level, the width of the central 
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FIG. 6. (Color Online) Upper panel: real part of the self en- 
ergy for the parameters U = — 2s f = 2, s c = 0.5, V 2 = 0.2, 
and W = 3.5, together with a piecewise linear fit (green 
dashed line). Middle panel: corresponding imaginary part 
of the self energy. Lower panel: corresponding local spectral 
function. 



Abrikosov-Suhl resonance is reduced, as is the kink en- 
ergy, cf. W = 4 in Fig. 

For oj > 0, the inflection point (hybridization gap) 
makes the identification of a kink more complicated. The 
most clear separation of kink (at oj ~ 0.004) and inflec- 
tion point (at uj ~ 0.01) is arguably obtained for W = 4 
in Fig. [3] But also for W = 6 two distinct features are 
discernible for oj > 0. 

Keeping the bandwidth fixed and modifying instead 
the hybridization strength, we show in Fig.Uthe real part 
of the self energy for the Bethc lattice. A kink for oj > 
is well visible, in particular for a smaller hybridization 
V, see right inset of Fig. 2] In the left inset of Fig. 2] we 
plot the kink frequency of the upper panel vs. the inverse 
hybridization strength 1/U 2 on a logarithmic scale. This 
reveals that the kink frequency shows the same exponen- 
tial dependence on 1/U 2 as the Kondo temperature. 

Let us note that compared to the Hubbard model for 
d-electron systems the kinks of the PAM are located at 
much lower energies (by about one order of magnitude). 
This stems from the very small quasiparticle renormal- 
ization Zfl *C 1) which is not surprising, since the inter- 
action strength is more enhanced in f-systems due to the 
confinement of electrons in the very narrow f-orbitals. 



Dispersion relation 

After finding kinks in the real part of the PAM self 
energy, let us now investigate if they have a similar 
influence on the effective dispersion relation as for 
the Hubbard model. Here, is defined as the max- 
imum of the fc-resolved spectrum A(h,oj) with respect 
to k, as in angular-resolved photo-emission experiments. 
Since the PAM is an effective two band model, it has 
two such dispersion relations for f and c(onduction) elec- 
trons. In Fig. [5] (a) and (b) the overall spectral density 
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FIG. 7. (Color Online) Specific heat of a Bethe-lattice PAM 
with U = — 2sf = 1, e c = 0.5, V 2 — 0.1 and various band- 
widths. 



Af(k,uj) + A c (k, uj) is plotted. The dispersion extracted 
from these spectral functions is depicted in Fig.[5](c) (red 
curve) At first glance, no kink feature is discernible. 

On the other hand, the kink in the real part of the self 
energy should directly reflect in a kink of the dispersion 
relation, whereas the imaginary part smears out the max- 
ima. This can be demonstrated by employing a linear fit 
to the real part of the self energy in Fig. [5] Taking into 
account the frequency range uj e [—0.05, 0.05], we obtain 
a Fermi liquid renormalization factor Zfl = 0.0689 for 
the slope around the Fermi level, and two renormaliza- 



tion factors Zgp 



0.202 and Z^ p = 0.312 for the slopes 



after the kink at u)*_ ss —0.005 and u>+ ~ 0.005, respec- 
tively; cf. Eq. ([2]). These renormalization factors are re- 
lated to the corresponding self energy slopes dT,/dui as 
Z = (1-dZ/du)- 1 . 

From the Fermi liquid renormalization Zfl and from 
Zq P the two dashed dispersion relations in Fig. [5] (c) are 
derived. These describe the NRG dispersion relations ac- 
curately around the Fermi level and for more negative fre- 
quencies. In-between there is a crossover from one curve 
to the other. This reflects the self-energy kink, which due 
to the already strong curvature of the non-interacting (or 
renormalized) dispersion does however not show up as an 
abrupt change of slope. 



Specific heat 

The kinks in the self energy can also be expected to re- 
flect as a change of the linear behavior of the specific heat. 
As in Ref. [1(3, we have employed the relation^ between 



low temperature entropy and spectral function, yielding 
the following conduction and f-electron contribution to 
the specific heat at temperature T: 



/oo 
dy 
-oo 



y 2 ev 



(ev + 1) : 



[A c (yT) 



A f(v T ) 
z(yT) 



— ReGf(yT)lmE(yT)] , 



(15) 



where Z(yT) = (1 — ReS(yT)) -1 is the renormalization 
factor for the /-electrons, and the summation over the 
two spin directions is accounted explicitly by the prefac- 
tor 2. The last term (second line) in Eq. (fT5j) also accounts 
for the imaginary part of the self energy and is beyond 
Ref. M 

Fig. [7] shows the specific heat for the Bethe lattice and 
different bandwidths (and density of states) of the con- 
duction electrons calculated according to Eq. (fT5"j) . At low 
temperatures, there is a linear increase of the specific heat 
as is to be expected for a Fermi liquid. The next domi- 
nant feature is a minimum found, e.g., at T 0.0006 for 
W = 3. The origin of this minimum is the hybridization 
gap which leads to a reduced number of states in the cor- 
responding energy interval. These two dominant features 
hide a more delicate kink feature which according to Eq. 
(fl5| should be present (note Z(yT) strongly depends on 
temperature), but is not discernible. 



III. CONCLUSION 

We have identified kinks in the real part of the self en- 
ergy of the periodic Anderson model, the arguably sim- 
plest model for f-electrons systems. The hybridization 
gap leads to an additional feature, in our case at uj > 0, 
making the clear identification of the kink more difficult 
than for the Hubbard model. The kink frequency fol- 
lows the same exponential dependence on the hybridiza- 
tion strength as the Kondo temperature. In contrast to 
the Hubbard model, the kink is difficult to detect visi- 
bly in the energy-momentum dispersion since the non- 
interacting energy-momentum dispersion has a strong 
curvature around the kink energy. Similarly, also in the 
specific heat, the fingerprint of the self-energy kink is less 
obvious because a stronger feature, associated with the 
hybridization gap, is superimposed. 
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